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Abstract 

In this paper we theoretically study electromagnetic reflection, transmission, and scattering prop- 
erties of periodic and random arrays of particles which exhibit both electric-mode and magnetic- 
mode resonances. We compare the properties of regular and random grids and explain recently ob- 
served dramatic differences in resonance broadening in the electric and magnetic modes of random 
arrays. We show that randomness in the particle positioning influences equally on the scattering 
loss from both electric and magnetic dipoles, however, the observed resonance broadening can be 
very different depending on the absorption level in different modes as well as on the average electri- 
cal distance between the particles. The theory is illustrated by an example of a planar metasurface 
composed of cut-wire pairs. We show that in this particular case at the magnetic resonance the 
array response is almost not affected by positioning randomness due to lower frequency and higher 
absorption losses in that mode. The developed model allows predictions of behavior of random 
grids based on the knowledge of polarizabilities of single inclusions. 

PACS numbers: 78.20.Ci, 42.25.Gy, 73.20.Mf, 78.67.Bf 
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I. INTRODUCTION 



Metamaterials are artificial composite materials which possess unusual electromagnetic 
properties not normally found in natural materials. Electromagnetic properties of nanos- 
tructured metamaterials in the optical range are one of the foci of interest in modern elec- 
tromagnetics. Traditionally, metamaterials and metasurfaces composed of small individual 
resonant inclusions are realized as periodical arrays. However, most recently, random or 
amorphous metamaterials start to attract attention, see Refs. 1 - 8 . This is due to novel tech- 
nological possibilities to manufacture amorphous structures cheaply and on a large scale, 
using advanced self-assembly techniques. In addition, effects of strong spatial dispersion 
(often undesirable) can be in some cases suppressed in disordered structures. It is generally 
accepted that the electromagnetic properties of both regular and random arrays of scatter- 
ers are quite similar if the distances between inclusions are electrically small. The main 
difference in electromagnetic response comes from scattering on the lattice inhomogeneities. 
This apparently results in additional loss in amorphous metamaterials, and for this reason 
regular metamaterial lattices have been the preferred choice if low-loss response is desired. 

However, it appears that in metamaterial structures exhibiting resonant responses in 
several modes, the effects due to position randomness of inclusions are more complicated. 
In a recent paper 1 by Helgert et al. reflection and transmission properties of regular and 
random (amorphous) planar arrays of cut-wire particles were studied both numerically and 
experimentally. Specially introduced position disorder of individual scatterers allowed to 
study the effect of distortion of periodicity on the electromagnetic response of the array. It 
was found that position randomness drastically affects the electromagnetic behavior at the 
electric resonance, but makes little impact at the array properties near the magnetic reso- 
nance of the particles. These results were validated by numerical simulations and confirmed 
in posterior work 2 . 

The authors of paper 1 put forward a hypothesis that the discovered dramatic difference 
between scattering properties in electric and magnetic modes is caused by difference in elec- 
tromagnetic interactions between particles in different modes. It was based on an observation 
that magnetic dipoles as well as electric quadrupoles do not generate tangential electric fields 
in the array plane, and it was assumed that this means that magnetic scatterers are not 
interacting with each other, so that the exciting field acting on a single particle is solely 
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the external illumination. On the other hand, the electric-dipole scatterers interact strongly 
and the exciting field is affected by positional disorder, which leads to resonance broadening 
and damping. However, from the duality principle it is known that in fact magnetic dipole 
particles interact via their magnetic fields exactly as strongly as electric dipoles interact via 
their electric fields, which means that the phenomenon discovered in Ref. 1 must have some 
other physical reasons. 

The goal of this paper is to study the phenomenon of resonance damping and broadening 
theoretically and explain the strong differences in resonance broadening in different resonant 
modes. To this end, we analytically study the effect of positional randomness on electromag- 
netic behavior of grids of resonant particles which can exhibit both electric and magnetic 
resonant responses. We introduce a simple model, which allows us to analyze the reflective, 
transmitting, and absorptive properties of multi-resonant grids, both in the regular and 
amorphous states. The theory is confirmed by numerical simulations using an example of 
the same metasurface as that studied in Ref. 1 . The results reveal the mechanisms of reso- 
nance broadening and damping in amorphous structures and explain the earlier discovered 
differences in the cases of electric (symmetric) and magnetic (anti-symmetric) resonances. 
Understanding physical phenomena which define the differences between effective electro- 
magnetic responses of regular and disordered metamaterials is urgently needed before the 
emerging amorphous metamaterials can find applications. Developing analytical models of 
disordered structures will allow the design and optimization of future composite materials 
with desired performance. 



II. ANALYTICAL THEORY OF PLANAR ARRAYS WITH ELECTRICALLY 
AND MAGNETICALLY RESONANT INCLUSIONS 

Let us consider an optically dense planar array of optically small resonant particles ex- 
cited by normally incident plane waves. We assume that the distance between the particles 
in the grid a is smaller than the wavelength. We are interested in the case when each 
particle exhibits both electric and magnetic responses, that is, both electric and magnetic 
moments are induced by local electric and magnetic fields, respectively. We also assume 
that bi-anisotropic magnetoelectric coupling is either forbidden due to the particle symme- 
try or it is negligible. Many widely-studied infra-red and optical metamaterial structures 
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like the cut- wire pairs considered in Ref. 1 belong to this class. In this paper we consider 
only electric and magnetic dipole moments of particles, neglecting quadrupoles and higher- 
order moments, concentrating on the influence of array randomness on the reflection and 
transmission coefficients. Relative strengths of dipolar and higher-order effects in cut-wire 
pairs have been analyzed in Ref. 9 . 

Assuming for simplicity that no cross-polarized dipole moments in the array plane are 
induced (the particles have the form of discs or squares, for example) and considering the 
excitation by normally incident plane waves, we can write the relations between the induced 
electric dipole moment p, magnetic dipole moment m, and the incident fields E inc and H inc 
as scalar relations 

p = a ee (E inc + (3eeP), m = a mm (H inc + (3 mm m) (1) 

Here a ee and a mm are the electric and magnetic polarizabilities of individual inclusions, 
respectively. Parameters (3 ee and (3 mm are called interaction constants and they measure 
contributions of the fields created by all other particles of the array into the local field 
E\ oc = E inc + f3 ee p exciting each particle (see e.g. Ref. 10 ). The interaction constants for 
electric and magnetic dipoles are related simply as 

ftmm o Pee (2) 

m 

where r] = y^tto/eo is the wave impedance of the surrounding space. Fields created by 
magnetic dipoles do not contribute to the electric local field exciting electric dipoles because 
the tangential component of the electric field of the magnetic dipole grid equals zero in 
the array plane. Likewise, fields scattered by electric dipoles do not excite magnetically 
polarizable particles positioned in the same plane. Most often, both moments are actually 
induced in the same particles, but the two modes have resonances at different frequencies. 

Next, we calculate the plane-wave electric fields created by the surface averaged electric 
current sheet J e = —^r and the magnetic current sheet J m = — %Ji ^r (the harmonic time 
dependence assumption is of the form e~ twt ): 

^ref = - 7T^e, H™ { = -^—Jm (3) 

E^ = - Vo H^, E Te f = E° cl + £™ f (4) 
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Here E^ e{ and E™ { are reflected electric fields created by the induced electric and magnetic 
currents J e and J m , respectively, and H™ f is the reflected magnetic field created by the 
induced magnetic current J m . Solving (1) for the induced dipole moments in terms of the 
incident fields and using (3) and (4) we find the reflection and transmission coefficients in 
the simple form 

^-T^-^ + -Km-Tl — I a o 2~ a ^ 

Emc 2a 2 - /3 ee 27] a 2 -i p mm 

"ee "mm 

T = 1 + R e - R m (6) 

Here we have used the plane-wave relation between the electric and magnetic incident fields: 
H inc = E inc /i] . The two partial reflections coefficients R e and R m correspond to the fields 
created by the induced electric and magnetic currents, respectively. Since (3 ee has the di- 
mension of l/(e a 3 ) and (3 mm has the dimension of l/(/zoa 3 ), it is convenient to multiply and 
divide the reflection coefficients by e a 3 or /i a 3 . The result is 

Re = ^ a' (7) 

2 e ° — R 
ik a 1 

"mm 

where ko = u^e fi is the wave number in the surrounding space. The normalized dimen- 
sionless interaction constants are the same for both electric and magnetic particles, and we 
denote them as (3: 

(3 = e a 3 /3 ee = /i a 3 /3 mm (9) 

Let us assume a simple Lorentz-type resonant response model of individual particles. This 
type of resonant response is very common and approximates very well the particle response 
near their resonances. Let us write down the inverse values of the normalized polarizabilities 
to make it easy to discuss the radiation loss factor: 

a ee V w oe - - iujT e J 6tt 

/loa 3 / A m \ ./end 3 . . 

' ( n ) 



Here T e ^ m model the dissipation losses in the particle (in respective modes), while the last 
imaginary term is due to the scattering (re-radiation of power) loss 10 . In case of regular or 



"totally random" (on the wavelength scale) grids there is no scattering loss, when the array 
period is smaller than the wavelength. In this case spherical-wave scattering from individual 
particles is suppressed by interactions between the particles in the array. Correspondingly, 
the imaginary parts of the interaction constants (3 ee and (3 mm contain terms proportional to 
fcg which compensate the corresponding terms in the inverse polarizabilities (see e.g. Ref. 10 ): 

Regular = MP) ~ ^ + ^ W 

The other imaginary term corresponds to the plane waves created by the surface-averaged 
currents. In case of amorphous (on the wavelength scale) arrays particles scatter individually, 
and there is no corresponding term in the interaction constants: 

/3amorph = Re(/3)+^ (13) 

In the quasi-static limit Re(/3) ~ 0.36 (see Ref. 10 ). 

Next, we substitute these interaction constants and the Lorentz particle polarizabilities 
(10) and (11) in the general formulas for the reflection coefficients (7) and (8). For regular 
or totally random (on the wavelength scale) arrays we get 

P _ -ha A , . 

eregular * 2 col -u 2 - <wr e - i k -fA e [ } 

_ k Q a Am 

2 u£ rn -id 2 - iujY m - i^A m 
Here Coq denotes the resonant frequency shifted due to interactions between the particles 
in the grid. In the quasi-static approximation for the real part of the interaction constant 
&oe,m ~ ^cL,™ ~~ 0.36A em . For amorphous grids we get 

_ - kpCL Ae , , 

^eamorph-? ^ . fc 3 3 , 

1 ul e -uj 2 - iujY e - -^A e - i<!fA e 
_ . k a Am , _s 

J*mamorph t n „i.3„3 , \*- ' ) 

Let us consider the case when electric and magnetic resonances occur at different frequen- 
cies. Then in the vicinity of one of the resonances the non-resonant moment varies weakly 
with the frequency and we can find a simple estimation of the resonant curve width (on the 
field-strength scale): 

9A -F 1 ^ m fio\ 

^^^e,m regular J- e,m ~r ~ _ V "/ 

fi 



for regular grids and 



^^e,mamorph J- e,m Z "I" ~ Z V J 

07T ^>0e,m * ^>0e,m 

for amorphous grids. 

We now see that if the condition 

^e,Jf^ + ^ » ^ (20) 

is satisfied, near the corresponding resonant frequency uo e ,m the effect of inclusion position 
randomness is negligible, and the response of regular and amorphous structures is nearly 
the same. Physically, this condition means that absorption (the first member of the left- 
hand side) and coherent plane-wave reflection (the second member on the left) dominate 
over scattering (the right-hand side term). The above relation shows that this is the case 
of high dissipative losses, low resonance strength, and small electrical size of the unit cell. 
Note that for the case of negligible absorption, this condition simply tells that scattering 
loss is negligible in random arrays if the distance between particles is optically very small 
(k 2 a 2 < 3tt). 

From the above results we can conclude that the effect of strong widening of the resonant 
curve of the electric-dipole mode and hardly any effect of array randomness on the magnetic 
mode discovered in Ref. 1 can be due to two reasons: 

1. At the frequency of the magnetic resonance the grid is practically homogeneous on the 
wavelength scale ("totally random"). Then the scattering term cancels out just like 
for periodical grids, and there is no difference in the resonant curve widths for regular 
and amorphous layers. 

2. At the magnetic resonance the particles are considerably more lossy and weaker excited 
than at the electric resonance, that is, (20) is satisfied near the magnetic resonance 
but not satisfied near the electric-mode resonance. 

III. EXAMPLE: ARRAYS OF CUT- WIRE PAIRS 

As an example we consider the cut- wire pair structure which was studied in Refs. 1,11,12 . A 
unit cell of the infinite regular array is depicted in Fig. 1. The dimensions are the same as in 




FIG. 1. (Color online) Geometry of the unit cell of the cut- wire array. 

Ref. 1 . The square lattice has the period of a = 512 nm along two transverse directions. The 
width of the cut-wire pairs in both lateral directions is W c = 180 nm. The height of each gold 
pair is H c = 30 nm. The gap between the two elements in each pair equals g = 45 nm and 
it is filled with a material with the relative permittivity equal to e r \ = 1.72. The structure 
is placed on top of a substrate with the permittivity of e r 2 = 1.5. The permittivity of gold 
is taken from Ref. 13 . 
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FIG. 2. (a): Electric polar izability of a single cut- wire pair 
and (b): Magnetic polarizability of the same particle 
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FIG. 3. (a): Electric polarizability of a single cut-wire pair in free space and (b): Magnetic 
polarizability of the same particle 

First we calculate the reflection and transmission coefficients for the regular array using 
the full-wave numerical simulator Ansoft HFSS. Substituting the numerical data for reflec- 
tion and transmission coefficients in (7) and (8) and using the quasi-static approximation 
for the real part of the interaction constants (3 ~ 0.36, we extract the polarizabilities of 
individual inclusions. The results are shown in Fig. 2, and it is apparent that the array has 
electric and magnetic resonances in different frequency regions, as expected. In fact the array 
is weakly bi-anisotropic due to the presence of the substrate (omega-type magnetoelectric 
coupling 14,15 ), which has been neglected in the theory and in the parameter extraction. We 
have checked that this approximation is valid by repeating the simulations and parameter 
extraction for the same array in free space. The results are presented in Fig. 3 and they 
show that this simplifying assumption is reasonable: The substrate effect is quite small. 
Behavior of the extracted electric and magnetic polarizabilities is very close to the canonical 
Lorentz-type resonant response. 

Scattering losses which appear in transition from regular to amorphous grids we model 
by the randomness parameter < r n < 1, where unity corresponds to the case where the 
scattering loss is completely compensated (regular array) and r n = means that the scat- 
tering loss is not compensated at all (amorphous array, each inclusion scatters individually). 
Transition from regular to amorphous state we model modifying the interaction constant 
(12) as follows: 

/3 = Re{/3}-r n ^ + ^ (21) 
which corresponds to a continuous transition from (12) to (13) with r n changing from unity 
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FIG. 4. (Color online) (a) - Amplitude of the reflection coefficient; (b) - Amplitude of the trans- 
mission coefficient; (c) - Phase of the reflection coefficient; and (d) - Phase of the transmission 
coefficient for grids with different randomness levels(r n ) 



to zero. It should be noted that for simplicity the randomness factor r n is assumed to be the 
same for both electric and magnetic interaction constants. Due to differences in resonant 
frequencies, this means that the same value of r n may correspond to somewhat different 
degrees of geometrical randomness of particle positions. 

Next we investigate how the reflection, transmission, and extinction change in transi- 
tion from regular to amorphous states, using the analytical formulas (7) and (8) with the 
extracted values of the polarizabilities and the interaction constant (21). Figs. 4 and 5 
show the randomness effects. One can see that the developed simple model gives very good 
agreement with the experimental and numerical data from Ref. 1 . Electrical response of the 
grid is strongly influenced by randomness, while close to the magnetic resonance there is 
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FIG. 5. (Color online) Absorption in the array in transition from regular to amorphous states 

almost no dependence on randomness. The reason for this phenomenon is the difference in 
the ratio of the absorption and scattering losses. Fig. 5 shows that in the periodical case 
the absorption is much stronger at the magnetic resonance than at the electric one. Higher 
losses are mainly due to larger imaginary part of gold permittivity, which is more than two 
times higher at the magnetic resonance: Im(e)io22nm ~ 3.2, Im(e)797 nm ~ 1.5. 

In addition, for the case of the grid in free space (Fig. 3) we have fitted the numerically 
extracted polarizability curves to the Lorentz model (10) and (11) and extracted parameters 
T em and A em . This allowed us to find the values in inequality (20). At the resonant 
frequency of the electric polarizability we find that the left-hand side equals 0.4+2.3 while the 
right-hand side equals 5.6. Scattering effects are clearly dominating and position randomness 
changes the array response quite significantly. At the resonant frequency of the magnetic 
polarizability the left-hand side reads 2.3 + 2, while the right-hand side equals 3. In this 
case the terms are of the same order and the randomness effect is much weaker. Note that 
condition (20) is a simple approximation which assumes that the two resonances are sharp 
and well separated. In this particular example, in the frequency region of the magnetic 
resonance the electric dipoles in fact give a significant contribution to the total absorption 
and coherent reflection. 

IV. CONCLUSIONS 

In this paper we have developed a simple model which explains the electromagnetic effects 
in transition from regular to random states of resonant particle arrays. We have derived 
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a general condition under which randomizing particle positions gives only negligible effects 
on the reflection and transmission coefficients and explained the earlier discovered dramatic 
differences in resonance damping for electric and magnetic modes of particles. We have 
also shown that the physical phenomena leading to the resonance damping in amorphous 
structures are the same for electrically or magnetically polarizable particles. The widening 
of the resonances takes place due to additional scattering losses, which are compensated in 
the case of electrically dense periodical grids. 

Studying transition to the amorphous state for a particular example of cut-wire pairs we 
have found that the reason for the much weaker resonance widening and damping in the 
magnetic mode is strong absorption in that frequency range. When scattering losses are 
much smaller than the dissipation losses, they make little impact on the total extinction. 
On the contrary, at the higher-frequency electric resonance scattering losses are stronger 
than the dissipation ones, which leads to strong resonance damping and distortion in the 
random case. In other situations, different transition effects in different resonant modes 
can also be caused by differences in the electrical size of the unit cell. In the considered 
example, the array period is comparable with the wavelength, thus, even for geometrically 
random positions of the particles with respect to the cell centers, the array cannot be made 
homogeneous on the wavelength scale. 

Our findings can have important implications in understanding the physical differences in 
electromagnetic responses of regular and amorphous structures, in design of various meta- 
material structures for such applications as subwavelength imaging, control of thermal radi- 
ation, microwave, terahertz and optical absorbers, and others. Using the developed model 
it is possible to predict and engineer the effects of randomness, relaxing conventional re- 
quirements on strong periodicity and make use of inexpensive self-assembly techniques in 
production of metamaterials. 
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